The nucleon axial charge in full lattice QCD 
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The nucleon axial charge is calculated as a function of the pion mass in full QCD. Using domain 
wall valence quarks and improved staggered sea quarks, we present the first calculation with pion 
masses as light as 354 MeV and volumes as large as (3.5 fm)^. We show that finite volume effects are 
small for our volumes and that a constrained fit based on finite volume chiral perturbation theory 
agrees with experiment within 7% statistical errors. 



INTRODUCTION 

The nucleon axial charge, f/^i, is a fundamental prop- 
erty of the nucleon, governing (3 decay and providing 
a quantitative measure of spontaneous chiral symme- 
try breaking in low energy hadronic physics. The ax- 
ial vector form factor, gA{<f'), is defined by the momen- 
tum space matrix element of the isovector axial current, 
^/x = qitj.lbfq: 

{N{p + q,s')\A^{q)\N(p,s)) - 

M(p-fg,s')^ [5^(9^)7^75 + 9p[q^)q^,l5\u{p,s). (1) 

The axial charge is the form factor at zero momentum 
transfer and is accurately measured from neutron (3 de- 
cay, gA = <7a(0) = 1.2695 ± 0.0029 (in units of the vector 
charge, gv)- It is also useful to write the axial charge in 
terms of the zeroth moments with respect to the momen- 
tum fraction x of the light cone quark helicity distribu- 
tions, Ag(a;) — q^{x) — qi{x). For each quark flavor g, 
these moments are defined and related to forward nucleon 
matrix elements of twist-2 operators as 

(1)a9 = / da;[Ag(j;) + Ag(a;)] 
Jo 

2MNS^{l)^q - {N{p,s)\q-i^-i^q\N{p,s)), (2) 

and the axial charge is given by gA = (1)au — (l)Arf- 

In addition to governing neutron fi decay, the axial 
charge is physically important for several other reasons. 
The Adler-Weisberger sum rulejj,i2i] relates the excess of 
g\ beyond the trivial value of 1 for a point, structureless 
nucleon to an integral over the difference between the ■n'^p 
and TT^p scattering cross sections, which is dominated 
by coupling to the A(1230) resonance. The Goldberger- 
Treiman relation(3|, gA — fTrd-wNN /Mn , relates it to the 



pion decay constant, pion nucleon coupling constant, and 
nucleon mass to an accuracy of five percent. Finally, the 
matrix elements in Eq. (jSJ directly measure the contri- 
bution of each flavor to the spin of the nucleon, and the 
flavor singlet combination S = (l)Au+(l)Ad+(l)As spec- 
ifies the fraction of the total nucleon spin arising from the 
spin of the quarks. Thus, a first principles calculation of 
the axial charge is an essential step in understanding key 
issues in hadron structure ranging from chiral dynamics 
to the spin content of the nucleon. 

The axial charge is the ideal starting point in the quest 
for precision lattice calculation of hadron structure for 
several reasons. It is accurately measured experimen- 
tally and the isovector combination (1)aii — (l)A(i has 
no contributions from disconnected diagrams, which are 
much more computationally demanding than the con- 
nected diagrams considered in this work. The functional 
dependence on both and volume is known at small 
masses from chiral perturbation theory (xPT) |3, ^ and 
renormalization of the lattice axial vector current can be 
performed accurately nonperturbatively using the five di- 
mensional conserved current for domain wall fermions. 
Thus, conceptually, it is a "gold plated" test of our abil- 
ity to calculate hadron observables from first principles 
on the lattice. In addition, since it is known to be partic- 
ularly sensitive to finite lattice volume effects that reduce 
the contributions of the pion cloud[^,0|, it is also a strin- 
gent test of our control of finite volume artifacts. 

LATTICE CALCULATION IN THE CHIRAL 
REGIME 

The computational challenge of precision lattice QCD 
lies in simultaneously approaching the continuum limit 
of small lattice spacing, the chiral limit of low quark 
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TABLE I: The number of configurations, pion mass, renor- 
malization factor, and axial charge for each of the measure- 
ments on two lattice volumes. 





20'^ X 32 lattice V 


= (2.5 fm)'^ 




# configs. 


(MeV) 


Za 


gA 


425 


76i (2) 


l.i296 (6) 


1.167 (11) 


350 


693 (3) 


1.1197 (6) 


1.153 (16) 


564 


594 (2) 


1.1085 (5) 


1.193 (16) 


486 


498 (2) 


1.0994 (4) 


1.173 (29) 


655 


354 (3) 


1.0847 (6) 


1.244 (58) 




28^ X 32 lattice V 


= (3.5fm)» 




# configs. 


(MeV) 


Za 


9A 


270 


353 (1) 


1.0840 (5) 


1.212 (59) 



masses, and the large volume limit. For nucleon struc- 
ture, where the pion cloud extends to large distances and 
plays a major role, the latter two requirements are par- 
ticularly severe. Instead of referring to the unobservable 
bare quark mass, since niq oc mj, it is convenient to ex- 
press the quark mass dependence of observables by their 
dependence on mj. To include the essential physics of 
the pion cloud, the box size must be large compared to 
the pion Compton wavelength, m~^, and the ultimate 
computational cost of a full QCD lattice calculation in 
such a box including dynamical sea fermions has mass 
dependence approximately m~^. 

To address this exceedingly demanding regime of light 
pion mass and large volume with current computational 
resources, we have used a hybrid combination of com- 
putationally economical staggered sea quark configura- 
tions generated using the so-called Asqtad improved ac- 
tion by the MILC collaboration!^ Q with domain wall 
valence quarks that have chiral symmetry on the lattice. 
Improved staggered quarks have been strikingly success- 
ful in precision calculations 10] and even predictions of 
heavy quark systems jlTj. Although there is not yet a 
proof that the continuum limit of the fourth root of the 
staggered determinant has the correct continuum limit, 
recent renormalization group arguments suggest that it 
does|l2j. Assuming that the staggered continuum limit 
is correct, the continuum limit of the hybrid calculation 
will also be correct. Both the Asqtad and domain wall 
actions have leading errors of order in the lattice spac- 
ing a. 

As tabulated in Tabled calculations were performed 
at five pion masses down to 354 MeV on 20^ x 32 lat- 
tices with a spatial volume V = (2.5 fm)^ and the light- 
est mass, which is most sensitive to finite volume ef- 
fects, was also calculated on a 28"^ x 32 lattice with a 
spatial volume V = (3.5fm)'^. The scale is set by the 
lattice spacing a — 0.12406 fm determined from heavy 
quark spectroscopy 13] with an uncertainty of 2%. To 
reduce the effect of dislocations at this lattice spacing 
and thus improve the chiral properties of the domain wall 



fermions, the fermion action was defined to include HYP 
smearing ^3 of the gauge fields, in which each link is re- 
placed by a sum of links within the unit hypercube. This 
HYP smearing changes the action by an irrelevant op- 
erator, is restricted in range to a single lattice spacing, 
and has coefficients defined to minimize fiuctuations of 
Wilson loops. 

The domain wall fermion [isl 1 16l | propagators are calcu- 
lated on a 5-dimensional lattice with fifth dimension L5 
and 5-dimensional mass M5. The physical quark fields, 
q{x, i), reside on the 4-dimensional boundaries with fifth 
coordinate 1 and L5 and have bare quark parameter 
{am)^^^ . The value M5 = 1.7 was chosen on the ba- 
sis of spectral flow analyses to optimize the evaluation 
of domain wall propagators, and L5 = 16 was chosen [l7j 
to ensure that the residual mass characterizing residual 
chiral symmetry breaking is always less than 10% of the 
physical quark mass for the pion masses in Table The 
values of {am)^^^ were tuned to reproduce the lightest 
pion mass obtained with the Asqtad action^^ and as 
shown in Ref . |0| , the resulting nucleon masses were also 
consistent with the Asqtad masses. 

Lattice matrix elements of the axial current, 
(^197^75 91-^)'°"''^'^; calculated from the ratio of three 
point functions to two point functions using the same 
smeared nucleon sources and methodology as in refer- 
ence The renormalization factor, Za, for the four 
dimensional axial current operator = qj^j^q is cal- 
culated using the five dimensional conserved axial cur- 
rent for domain wall fermions A^, by the relation[l^ 
{A^{t)A,m ^ ZA{A^{t)A,m. 

LATTICE RESULTS 

The lattice renormalization constants, Za, and 
values of the renormalized axial charge, gA ~ 
ZA{N\qjf,-f5q\Ny''**''=<', are tabulated in Tabled for the 
five pion masses and two volumes considered in this work. 
Measurement of gA to 5 % accuracy at 171-^= 354 MeV 
is a major achievement of this work, requiring 655 con- 
figurations on the 20'^ lattice and correspondingly fewer 
on the larger lattice where volume averaging improves 
statistics. 

The dependence of gA on the pion mass HH Sim, 
and on both the mass and volume 0, is known from 
xPT, where at sufficiently low pion mass, the functional 
form follows from the symmetries of QCD and the coeffi- 
cients are low energy constants of QCD that in principle 
can be calculated on the lattice. We have used the xPT 
form of Ref. Q and verified that it agrees closely numer- 
ically with the other corresponding infinite volume [5ll2(]| 
and finite volumeQ results for our masses and volumes. 
Our xPT fit includes six constants: (the pion decay 
constant), mA —rriN (the Delta - nucleon mass splitting), 
9 A, 9 A A, and gNA (the couplings of the axial current be- 



T ' 1 ' 1 ' 1 ' 1 ' 1 , . r ' 1 ' 1 ' r 




1 I 1 I 1 I 1 I 1 I 1 I 1 I n\ 1 I 1 I 1 I 1 L_ 

0.1 0.2 0.3 0.4 0.5 0.6 0.2 0.4 0.6 0.8 

m^' (Gev') (Gev') 



FIG. 1: Nucleon axial charge qa as a function of the pion 
mass. Lattice data are denoted by squares (smaUer volume) 
and a triangle (larger volume), the lowest smaller volume 
point is displaced slightly to the right for clarity, and ex- 
periment is denoted by the circle. The heavy solid line and 
shaded error band show the xPT fit to the finite volume data 
evaluated in the infinite volume limit, and the lines below it 
show the behavior of this chiral fit in boxes of finite volume 
L^, as L is reduced to 3.5, 2.5, and 1.6 fm respectively. 



tween two nucleons, two Deltas, or a nucleon and Delta, 
respectively), and a counterterm C. 

Figure n shows the lattice data and our fit to it using 
finite volume xPT. The function was fit to each 

data point at the corresponding mass and finite volume, 
and the parameters of the fit were then used to determine 
the infinite volume axial charge. In the absence of lattice 
calculations of gA at still lower pion masses, it is not 
presently possible to do a complete extrapolation from 
lattice measurements alone. Hence, following Ref. |20| . 
we performed a constrained fit and the heavy solid curve 
is determined by setting /tt, ttia — tojv, and gjvA to their 
physical values '2(11 and performing a least squares fit for 
gAi 5AA, and C. The error band arising from this three 
parameter fit is shown in Fig. ^ and the resulting value 
for the axial charge at the physical pion mass is (?A(™7r — 
140 MeV) = 1.212 ±0.084. Given the smooth behavior of 
the chiral fit and the small magnitude of the chiral logs, 
it is clear that the extrapolation for the axial charge is 
quite benign, and that the lattice data is extrapolating 
convincingly towards experiment. 

Although there has been concern that qa is particu- 
larly sensitive to finite volume eS'ects0, 0, H^l, Fig. 
shows these efl[ects are well under control and introduce 
negligible errors for our volumes. The light curves show 
the behavior expected from xPT in volumes with L of 
3.5, 2.5, and 1.6 fm. Note that at the lightest mass, our 
3.5 and 2.5 fm results are statistically indistinguishable, 
consistent with the xPT change of less than 1 %, and 



FIG. 2: Comparison of all full QCD calculations of qa, as 
described in the text. The solid line and error band denote 
the infinite volume xPT fit of Fig. 1 and its continuation to 
higher masses is indicated by the dotted line. 



that the corrections applied in correcting our data from 
2.5 or 3.5 fm to infinity in the xPT fit are quite small. At 
heavier masses, although the truncated xPT expansion is 
not quantitatively reliable, the finite volume effects are 
physically suppressed. The order of magnitude of the 
corrections from 1.6 to 2.5 fm is also consistent with the 
fact that quenched calculations!^ have shown that in- 
creasing the box length from 1.2 fm to 2.4 fm increases 
gA by the order of 10% for pion masses ranging from 550 
to 870 MeV, and unquenched calculations for 770 MeV 
pionsV| have shown that increasing L from 1.1 to 2.2 fm 
increases gA by 20%. 

In addition to the statistical error arising from fitting 
the parameters gA, 5a A, and C, several systematic errors 
may be estimated. The three constrained parameters can 
be calculated directly on the lattice and the linear re- 
sponse of our chiral fit to varying each of of them shows 
very weak dependence. Calculation of on our lattices 
yields 92.4 MeV ± 3%, corresponding to an error in gA 
of 0.10%, and a rough calculation of toa — rriN, (which 
can be improved) with 18% uncertainty corresponds to 
an error of 0.29%. Although we have not yet calculated 
gNA, it should be calculable to 20%, corresponding to 
an error in gA of 0.13 %. Thus the total error from the 
constrained parameters is much less than a percent. The 
error in the lattice scale, which can shift all masses by 2%, 
will induce a negligible effect since the curve in Fig.^is so 
flat. An alternative lattice renormalzation method based 
on calculating the ratio of the axial and vector charges, 
which should have the same renormalization constant in 
the chiral limit, yields discrepancies less than 2% for the 
heaviest masses and statistically indistinguishable results 
at lighter masses, suggesting that the renormalization er- 
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ror is less than a few percent. 

Finally, we compare our results with the world's sup- 
ply of unquenched lattice calculations of in Fig. |21 
The dotted line shows the continuation of our xPT fit 
to higher masses to guide the eye, but is not quanti- 
tatively reliable due to the limitations of xPT at such 
large pion masses. O ur p revious calculations |0 using 
SESAM configurations [22j with dynamical Wilson quarks 
and spatial box size 1.5 fm are shown for the heaviest 
three quark masses. These are consistent with the cur- 
rent hybrid calculations within statistics, although the 
lowest mass point may have a small systematic shift 
downwards due to the small volume. The three points 
by the RBCK collaboration^24iJ using dynamical domain 
wall fermions in a 1.9 fm box are also consistent with 
our present results within statistics. The initial results 
by the QCDSF/UKQCD collaboration^ using improved 
Wilson dynamical fermions in boxes ranging from 1.5 to 
2.2 fm are systematically lower than the results of the 
other three calculations, although finite volume effects 
and differences in renormalization away from the chiral 
limit may account for the discrepancy. 



CONCLUSIONS 

In summary, we have calculated gA hi full QCD in the 
chiral regime. The hybrid combination of improved stag- 
gered sea quarks and domain wall valence quarks enabled 
us to extend calculations to the lightest mass, 354 MeV, 
and largest box size, 3.5 fm, yet attained, and to ob- 
tain statistical accuracy of 5% with negligible error from 
volume dependence. Chiral perturbation theory implies 
mild dependence on the pion mass, and a three param- 
eter constrained fit yields an excellent fit to the data 
and generates an error band of size 7% at the physical 
pion mass which overlaps experiment. Thus, this calcu- 
lation represents a significant milestone in the quest to 
calculate hadron structure from first principles. Building 
on this result, future challenges will include extension of 
the range of pion masses to include 300 and 250 MeV 
to enable an unrestricted determination of the param- 
eters of xPT from lattice data, generalization to par- 
tially quenched calculations and the application of par- 
tially quenched hybrid xPTji^, and evaluation of dis- 
connected diagrams to enable calculation of the quark 
spin structure of the nucleon. 
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